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In a previous paper, we studied the ergodic properties of an Euler scheme of a 
stochastic differential equation with a Gaussian additive noise in order to approximate 
the stationary regime of such equation. We now consider the case of multiplicative 
noise when the Gaussian process is a fractional Brownian Motion with Hurst parameter 
p ^ H > 1/2 and obtain some (functional) convergences properties of some empirical 

^ measures of the Euler scheme to the stationary solutions of such SDEs. 
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1 Introduction 



Stochastic Differential Equations (SDEs) driven by a fractional Brownian motion (fBm) 
have been introduced to model random evolution phenomena whose noise has long range 
dependence properties. Indeed, beyond the historical motivations in hydrology and telecom- 
munication for the use of fBm (highlighted e.g in [21]), recent applications of dynamical 
systems driven by this process include challenging issues in Finance [11] , Biotechnology [23] 
^ ; or Biophysics [El [16]. As a consequence, SDEs driven by fBm have been widely studied 

in a finite-time horizon during the last decades, and the reader is referred to |23|. 0] for 
nice overviews on this topic. 

In a somehow different direction, the study of the long-time behavior (under some 
stability properties) for fractional SDEs has been developed by Hairer (see [12j [13]) and 
Hairer and Ohashi |14j . who built a way to define stationary solutions of these a priori non- 
Markov processes and to extend some of the tools of the Markovian theory to this setting. 
See also [U 02 [10] for another setting called Random dynamical systems. The current 
article fits into this global aim, and starts from the following observation: the knowledge 
of the stationary regime being important for applications and essentially inaccessible in an 
explicit form, we propose to build and to study a procedure for its approximation in the 
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case of SDEs driven by fBm with a Hust parameter H > 1/2. This paper is following a 
similar previous work for SDEs driven by more general noises but in the specific additive 
case (see [3]). 

More precisely, we deal with an Revalued process (X t )t>o solution to the following 
SDE 

dX t = b{X t )dt + a{X t )dB* (1) 

where b : R rf — > R rf and a : M. d — > M^g are (at least) continuous functions, and where 
M^o is the set of d x q real matrices. In ([1]), (B^)t>o is a q-dimensional .ff-fBm and for 
the sake of simplicity we assume ^ < H < 1, which allows in particular to invoke Young 
integration techniques in order to define stochastic integrals with respect to B. Compared 
to [3] we handle here a fairly general diffusion coefficient a, instead of the constant one 
considered previously. Classically the noise is called multiplicative in this setting, whereas 
it is called additive when a is constant. 

Under some Holder regularity assumptions on the coefficients (see Section [2] for details), 
(strong) existence and uniqueness hold for the solution to ([!]) starting from xq G R d . Clas- 
sically for any stochastic differential equation, a natural question arises: if we assume that 
some Lyapunov assumptions hold on the drift term, does it imply that (Xt)t>o has some 
convergence properties (in distribution or in a pathwise sense) to a steady state when 
t -)• +oo ? 

This question implies in particular to define rigorously a concept of steady state. For equa- 
tion (pQ), this work has been done in p3]: using the fact that, owing to the Mandelbrot 
representation, the evolution of the fBm can be represented through a Feller transition on 
a functional space S, the authors show that a solution to (pQ) can be built as the first coor- 
dinate of an homogeneous Markov process on the product space M. d xS. As a consequence, 
stationary regimes associated with (JTJ) can be naturally defined as the first projection of 
invariant measures of this Markov process. Furthermore, the authors of [14] develop some 
specific theory on strong Feller and irreducibility properties to prove uniqueness of invari- 
ant measures in this context. 

In the current article, our aim is to propose a way to approximate numerically the station- 
ary solutions to equation ([1]). To this end, we study some empirical occupation measures 
corresponding to an Euler type approximation of ([I]) with step 7 > 0. We show that, 
under some Lyapunov assumptions, a (pathwise) convergence to the stationary solution 
of the discretized equation (denoted by v 1 ) holds and that, when 7 — > + , v" 1 converges in 
turn to the stationary solution of (P). This approach is the same as in [3]. However, the 
introduction of multiplicative noise has some important consequences on the techniques 
for proving the long-time stability of the Euler scheme. In particular, the main difficulty is 
to show that the long-time control of the dynamical system can be achieved independently 
of 7. In [3j, this problem has been solved with the help of explicit computations for an 
Ornstein-Uhlenbeck type process. Because the noise is multiplicative the computations 
of [3] are not feasible anymore and we use specific tools to obtain uniforms controls of 
discretized integrals with respect to the fBm. Before going more precisely to the heart of 
the matter, let us mention that the numerical approximation of the stationary regime by 
occupation measures of Euler schemes is a classical problem in a Markov setting including 
diffusions and Levy driven SDEs (see e.g. [281 E3 CE1 ESI E5J [26]). 

2 Framework and main results 

For some fixed H G (|, 1), we consider (f2, F, P) the canonical probability space associated 
with the fractional Brownian motion indexed by R with Hurst parameter H. That is, 
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0, = Cq(R) is the Banach space of continuous functions vanishing at equipped with the 
supremum norm, T is the Borel sigma-algebra and P is the unique probability measure on 
Q such that the canonical process B = {Bt = (Bj, . . . , Bf), t G R} is a fractional Brownian 
motion with Hurst parameter H. In this context, let us recall that B is a g-dimensional 
centered Gaussian process such that Bq = 0, whose coordinates are independent and 
satisfy 



E 



B\-B\' Z 



t-s\ 2H , for s,teR. (2) 



In particular it can be shown, by a standard application of Kolmogorov's criterion, that 
B admits a continuous version whose paths are Holder continuous for any 9 < H. 

Let us be more specific about the definition of Holder spaces of continuous functions. 
Namely, our driving process B lies into a space C e defined as follows: we denote by 
C 9 (M+,lR d ) the set of functions / : R + -> R d such that 

VT>0, ||/|| 0)T = sup l^-^l <+oo, 

0<s<t<T (t-sf 

where the Euclidean norm is denoted by | . |. We recall that C e (R+,R d ) can be made into 
a non-separable complete metric space, whenever endowed with the distance 5g defined by 



sup \\f(t)-g(t)\\ + \\f-g\\ eiN 

0<t<N 



where x A y = min(x, y) \/x, y G R. However, since separable spaces are crucial for conver- 
gence in law issues, we will work in fact with a smaller space C e (M+, R d ): we say that a 
function / in C e (M + ,R ci ) belongs to C e (M + ,R d ) if 

wT^n r* x\ !/(*) ~ f( s )\ S-+0 K n 
VT>0, uj0 t T{f,S):= sup — rg > 0. 

0<s<t<T,0<\t~s\<8 I* — s l 

C (M+,lR d ) is a closed separable subspace of C e (M + ,M (i ). 

We recall now some results on the existence and uniqueness of the solutions of the 
stochastic differential equation ([1]) starting from a deterministic point a. Let us suppose 
that b is Lipschitz continuous and that a is (1 + a)-Lipschitz with a > i — 1,. We recall 
that for < a < 1, a is (1 + a)-Lipschitz on R d if it is a C 1 function and if the following 
norm is finite: 

, ■ i , \ 1 1 \Da(x) — Da(y)\ , . 

\\a\\ 1+a = sup \\Da(x)\\+ sup ^ (3) 

\x-y\<8 \ x — y\ 



Then for any deterministic function B £ C e (R+,R q ) with 9 > ^, any xq G R d and any 
(1 + a)-Lipschitz map a, there exists a unique solution denoted by &(xq, B) in C 9 (R+, R d ) 
of 

dX t = b(X t )dt + o{X t )dB t , X = x . (4) 

See [U [20] for proofs and discussions of this kind of deterministic ordinary differential 
equations. Please note that by definition 

<f>(x ,B) t = xo+ I b{$(x ,B) s )ds+ [ a(<S>{x ,B) s )dB s , 
Jo Jo 

where the integrals are Riemman Stieljes integrals. Moreover the so-called Ito map $ is 
continuous from R d x C e (R + ,M ,? ) into C e (M + ,M 9 ) (see Proposition 47 in [1]). Actually 
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the solution of the stochastic differential equation (pQ) such that Xq = xq is nothing but 
$>(xq, B H (u)), where the fractional Brownian motion B H has been first sampled. 

Let us now introduce a long-time stability assumption (C) concerning the stochastic dif- 
ferential equation (flj where b : M. d — > M. d is a continuous function. 

Let £Q(M. d ) denote the set of Essentially Quadratic functions, that is C 2 -functions 
V : R d -> (0, oo) such that 

liminf > o, |VV| < CW and D 2 V is bounded. 

[a:[-++oo \Xf 

Note that any element V G £Q(M. d ) is continuous, and thus attains its positive minimum 
v > so that, for any A, r > 0, there exists a real constant C A r such that A+V r < C A r V r . 

(C): The map a is assumed to be a bounded Lipschitz continuous function. Moreover we 
suppose that there exists V G £Q(M. d ) such that 

(i) for a given C > we have \b(x)\ 2 < CV(x) Vx G M d . 

(ii) and for /3 G R and a > 0, the following relation holds: 

VxGK ', (VF(x),6(x))</3-qF(i). 

We now turn to the definition of a stationary solution of ([T|). 

Definition 1. For 1/2 < 9 < H < 1, a stationary solution i/ to flU is a probabil- 
ity measure on C e (M + , R d ), which is the first projection of a probability measure \i on 
x C 9 (M,1R 9 ), satisfying the following conditions. Let us denote by (X t ,B t ) the 
coordinate process on C e (M + ,lR d ) x M 9 ). The projection of fi on C e (R,R 9 ) is the law 

of a ff-fractional Brownian motion indexed by R. For any T > 0, we consider [i? the 
restriction of (i on C e ([0, T], R d ) x C e ([0,T],R q ). Then we assume that 

X = <J>(X , {Bf) t >o), Ht almost surely (5) 

and that the distribution of {Xt)t>o is strictly stationary under /i. A stationary solution is 
called adapted, if for < s < t the processes (X s )o< s <t and (B H (s)) s >t are conditionally 
independent given (B H (s)) s <f 

This definition is very similar to the definition in [T3] that comes from Stochastic Dynam- 
ical Systems. 

Let 7 be a positive number, we will now discretize equation ([1]) as follows: 

Y? = yj 7 + (t- n 7 )6(yj 7 ) + a(Y^)(B t H - B*) Vt G [n 7 , (n + 1) 7 ). (6) 

We set 

t 7 = max{jk,jk < t, k G N}. 

In fact, we will usually write i instead of t„ in the sequel. The discretization of ([TJ can 
also be introduced with the following discretization <3? 7 : M x C e (M+,M 9 ) (->• C e (]R + ,IR d ) of 
the ltd map : 

&{x , B) t := x + / 6(<*> 7 (x , B)s)ds + I a(^(x , B) )dB s . (7) 
Jo Jo 
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Please note that the definition of $ 7 does not involve any Riemann integration but only 
finite sums and that 

TH = &(Y Q \ (Bf ) t >o)f (8) 

Then, we call adapted solution of ([6]) a distribution v 1 on C e {R+,R d ) defined as in Def- 
inition Q] replacing (p>J) by ([8]) . We will say that v" 1 is stationary if it is invariant by the 
shift maps {6 kl ) km , where t : C {R + ,R d ) -»■ C e (M+,IR d ) is such that (0 t (w))„ = u t +s, 
for t > 0. Note that in this definition, there is a a slight abuse of language since we do 
not require the invariance by the shift maps Q% for every t > 0, but only when i = k'j, k G N. 

Let us introduce the following uniqueness assumption for v~' and v: 

(S 7 ) (7 > 0): There is at most one adapted stationary solution to (P) (resp. to ©) if 
7 = (resp. if 7 > 0). 

For (S°), we refer to Theorem 1.1. of [14J. When 7 > 0, we have the following proposition: 

Proposition 1. Let H G (1/2, 1). Assume that d = q and that b and a are C 2 -functions. 
Assume that a is invertible and that sup,j, gR d a~ l {x) < +00. Then, (S 7 ) holds for every 
7 > 0. 



The proof, which is an application of (13], is done in the appendix. 

Let us now focus on the construction of the approximation: we denote by (X t 7 )t>o the 

■7 




continuous-time Euler scheme defined by: Xl =i£l' ! and for every n > 



X? = X 7 + (t - n 7 )6(X 7 ) + a(X^)(B? - B*) Vt G [n 7 , (n + 1) 7 ). (9) 



The process (X^ )t>o is a solution to ([6]) such that Xq = x. In order to alleviate the 
notations and, when it is not confusing, we will usually write Xf instead of X^ . Now, we 
define a sequence of random probability measures (P(")(a;,da)) n >i on (7 e (R+,]R d ) with 
6 < H by 

1 n 

^ (n ' 7) (-^«) = -E^ Vl)+ H(^) 

fc=l 

where 5 denotes the Dirac measure and for every s > 0, := (Xj +t )t>o denotes the 

s-shifted process. 

For t > 0, the sequence (7^ (w, dy)) n >i of "marginal" empirical measures at time t on 
M d is defined by 

1 n 

n (n ' 7) (u;,d y ) = -^^ U+ J.)(dy). 

n A — ' 7(fc-i)+* v ' 

k=l 

We are now able to state the main theorem of this article: 

Theorem 1. Let 1/2 < 6 < H < 1. Assume (C) and (S 7 ) /or every 7 > 0. T/iera, 
(%) There exists 70 > such that for every 7 G (0,7o), 



V^'^^uj, da) ^ C ^ z; 7 (da) a.s. where n 



-00 



and where v 1 denotes the unique adapted stationary solution to ©. 

(ii) Moreover, if b is Lipschitz continuous, a is (1 + a)-Lipschitz with a > jj- — 1 and if 
(S°) holds, then, 

lim lim "P ( - n ' 7 ^(aj, da) = z^(da) a.s. 

7— >0 n— >+oo 

where the convergence is again for the weak topology induced by C e (R+, R d ) and v denotes 
the unique adapted stationary solution to dU). 
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Remark 1. Note that some extensions can be deduced from the proof of this theorem. 
First, when uniqueness fails for the stationary solutions, the preceding result is replaced 
by 

Theorem 2. 1. Assume (C). Then, there exists 70 > such that for every 7 £ (0,70), 
(P^ n ' 7 ^(w, da)) n >i isa.s. tight on ID(IR + ,IR rf ). Furthermore, every weak limit is a stationary 
adapted solution of (|6j). 
2. Assume (C) and set 

^°°' 7 H := {weak limits of (P (n ' 7) (w, da))}. 

Then, there exists 71 G (0,70) such that (W°' 7 (o;)) 7 < 7l is a.s. relatively compact for 
the uniform convergence topology on compact sets and any weak limit when 7—7-00/ 
(Z//°°' 7 (a;)) 7 < 7l is an adapted stationary solution of ([T]). 

Remark 2. From the very definition of weak convergence, the preceding assertions im- 
ply that the convergence of {T ><ynn \uj, da)) n ^ holds for bounded continous functionals 
F : C e (R+,R d ) -> R. In fact, this convergence can be extended to some non-bounded 
continuous functionals F : C e ([0, T], R d ) — > R. Actually, setting G(a) = sup ig [ 0i T] V(at), 
we easily deduce from inequality (fT2|) of Proposition [3] and Proposition |4] that 

sup limsup V^\lo,G p ) < +00 a.s. 

7<70 n— ¥+00 

for every p > 0. By an uniform integrability argument, it follows 

Proposition 2. The convergence properties of ' (T , ^ n '" / ^ (oj , da)) extend to continuous func- 
tionals F : C ([O,T],R d ) -> R such that for every a 6 C e ([0, T],R d ), 

\F(a u < t<T)\<C sup V p {a t ) 

te[o,T] 

with T > and p > 0. 

Remark 3. A third natural extension of Theorem [1] consists in handling the case of an 
irregular fractional Brownian motion B with Hurst index 1/4 < H < 1/2. This extension 
is presumably within the reach of our technology on differential systems driven by fBm, but 
requires a huge amount of technical elaboration. Indeed, to start with, equation ([1]) has to 
be defined thanks to rough paths techniques whenever H < 1/2, and we refer to [9] for a 
complete account on rough differential equations driven by Gaussian processes in general 
and fractional Brownian motion in particular. More importantly, as it will be observed in 
the next sections, our main result heavily relies on some thorough estimates performed on 
the discretized version (|6|) of equation ([1]). When H > 1/2 this discretization procedure is 
based on an Euler type scheme, but the case H < 1/2 involves the introduction of some 
Levy area correction terms of Milstein type (see [6j) or products of increments of B if 
one desires to deal with an implementable numerical scheme (cf. [7]). This new setting 
has tremendous effects on the proof of Propositions [3] and [U For sake of conciseness, we 
have thus decided to stick to the case H > 1/2, and defer the rough case to a subsequent 
publication. 

The sequel of the paper is built as follows. The three next sections are devoted to 
the proof of Theorem [TJ In Section [3l we prove some preliminary results for the long- 
time stability of (V^ n,1 \oj, da)) n , when 7 > 0. It is important to note that the controls 
established in this section are independent of 7 in order to obtain in the sequel a long- 
time control that does not explode when 7 — > 0. Then, in Section U we obtain some 
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tightness properties for (V^ n '^(uj, da)) (in n and 7) and, in Section [5l we prove that the 
weak limits of this sequence are adapted stationary solutions. Eventually, in Section [6l 
we test numerically our algorithm for the approximation of the invariant distribution of a 
particular fractional SDE. 

Note that in the proofs below, every non-explicit constant is denoted by C and may change 
from line to line. 



3 Evolution control of (X^) in a finite horizon 

The main aim of this part is to obtain a finite-time control of V(Xj) in terms of V(Xq) 
which is somewhat contracting and which is independent of 7. This is the purpose of the 
first part of Proposition [3] below. In order to obtain some functional convergence results, 
we state in the second part a result about the finite-time control of the Holder semi-norm 
of JT. 

Proposition 3. Let T > 0. Assume (C). Then, 

(i) For every p > 1, there exist 70 > 0, p G (0, 1) and a polynomial function P Pj q : R — >• R 
such that for every 7 G (0, 70], 

V>(XZ) < P V*(x) + P p , e (\\B H \\e, T ) + CQP(B?, < t < T), (10) 

where Q 7 is defined for every (w(t))te[o,T] by 

LfJ 

Qj((Ht))te[o,T]) = E \ w ( k ^ ~ w « k ~ ^\ 2 - ( H ) 
fc=i 

sup V(^7) < C(V p (x) + P p ,o(\\B H \\o tT ) + Q p 1 (Bi I ,0 < t < T)). (12) 
te[o,T] 

^ For ei;ery 6* G (|,-ff), T > 0, and 7 G (0,70] 

sup '^"f/ 1 < C T (F(x) + P fl (||B ff || fl)T ) + Q 7 (B? , < t < T)) 

0<s<i<T (t - S) 

where P is another real valued polynomial function. 

The proof of this result is achieved in Subsection l3.2i Before, we focus in Subsection 13. II 
on the control of increments of some discretized equations with non-bounded coefficients 
driven by B H . 

3.1 Technical Lemmas 

Let us recall that, for every t > 0, t^ = 7max{fc G N,jk < t}. In the sequel, we will 
usually write t instead of 1L. 

In the following lemmas, we will use the following notation: for any element (x(t))t>o 
of C(M+,R d ) and T > 0, 9 > 0, 7 > 0, we define 

,, u s,t \x{v)-x{u)\ 

IfL-v = su p 



~<U<V<t [V-y ~ U-y) 



s< 1 u< 1 v< 1 t V— 7 —7 

where we set by convention § = 0. 



Lemma 1. Assume that b is a sublinear function, i.e. that there exists C > such that 
for every x G \b{x)\ < C(l + |x|). Then, for every T > 0, there exists a constant C > 
such that for every s,t G [0, T] wii/i s <t, for every 7 > 0, for every 8 G (0,H) 



|*7| < (|X2I + + ||^||^(t-#J exp(C7(t-s))) 

w/iere 

Z? = fa{X])dB». 
Jo 

Proof. First, from the very definition of (*7)f>o> we have for every s,t G [0, T] with s < t: 

*7 = XI + f~b{Xl)du + Z 7 - Z2_ (13) 
The function 6 being sublinear, we deduce that 

1*71 = 1*71 + ||zi#C£ - i) e + c Ai + \xi\)du. 

J s_ 

Setting g s (v) = \Xg+v\i it follows that for every v G [0, t — s], 

/•I! 

5s(^) < a + C I g s {u)du 



with a = |*7l + ll^lle \{t — s) e + C{t — s). The result follows from Gronwall's lemma. □ 

The control of S^-integrals is usually based on the so-called sewing Lemma (see e.g. 
[Jl [8] ) which leads to a comparison of J s f(x u )dB^ with f(x t ){Bf - Bf). The following 
lemma can be viewed as a discretized version of such results: 

Lemma 2. Let 70 > and (fj)-ye(0,j o ] ^ e a f am ^y °f functions from M + x M. d to M^g 
such that there exists r > such that for every T > ; there exists Ct > such that 
V7 G (0,7o], 

V(s,x),(t,y) G [0,T]xR d , \\f^t,y)- f 7 (s,x)\\ < C T (l+\x\ r + \y\ r )(\ts\ + \y-x\). (14) 
Let (H^)t>o be defined by 

t 



vt>o, W t = [ / 7 (s,*7)dsf 

J 



Then, for every 9 G (^,H), for every T > 0, f/iere exists CV > such that for every 
7 G (0, 70], for every for every < s < t < T , 

|i?7 -hi- f y (s,x2)(B t H - b* )\ < c T {t_ - sf\{i + |*2r +1 + {\\z^py^)\\B H \\ e , T . 

Proof. Let s, t G [0, T] with s = ^i and i = 7J with £ < j. We build the sequences (t, )f =0 
by: Tq '' = s, = i and for every k > 1, for every / G {0, . . . , 2 fc }: 



r (*) 



Tj 11 / is even 



7 2 i(r& 1) +-£7 1) ) 

2 2 



(15) 

if / is odd. 



(k) 

When I is odd, r, is in fact a discretization point which is located in the interval 
[t^ , Till More precisely, it can also be defined as the largest discretization point on 
the left of the middle of this interval. It is important to remark that for k large enough, 
every discretization point between s and t is covered by one of the r/ fc) 's, I £ {0,...,2 fc }. 
Then, we set 

K = Kij := inf jjfe E N, {7m, i < m < j} C {rf , < I < 2 fc }} . 

Before going further, let us check that there exists a constant C > (C = 2 is a suitable 
choice) such that for every 7 > 0, for every s, i S [0, T] 

VA;e{0,K}, V/G{l,...,2 fc }, r« - rffi < C^. (16) 

If = 0, this point is clearly satisfied. Assume now that K > 1. Denoting by Rjj^ the 

number of discretization points (strictly) between 7], and we can check by induction 
that 

\fk G {0, . . . , if}, max{nJJ^. 1 , < Z < 2 fc - 1} - mm{njff +1 ,0 < I < 2 k - 1} < 1. 
In other words, 

Vfc G {0, . . . , K}, maxjrjj - r/ fc) , < / < 2 k - 1} - mm{T^\ - r/ fe) , < I < 2 fc - 1} < 7. 
Using that ^/ 2 =i ~~ T /-i = £ — s, we first deduce that 



.(*) J*) 



t - s 



Owing to the definition of if, the second consequence is that {r^ K )/ is (strictly) increas- 
ing and thus that min 1< ;< 2 if-i (r^ K ^ — r,^ ) > 7. This implies that 

, t-s 

and ()16p follows. 

Now, let us return to the proof of the lemma. Using that at step K, the discretization 

(K) 

steps are covered by the r, s, we obtain: 



Hi ~ H] = f(Tl K \x\ K) ){B» K) - B« K) , 

~ l=Q r, r l+1 



Thanks to (|15p . we have for every k £ {0, . . . , K - 1} and I £ {0, . . . , 2 k - 1}: 

B ( fc) - B {k) = (i?5fe+i) ~ B J.k+i)) + (Bjk+i) ~ Bj k+1) ). 
'1+1 '1 '21+2 '21+1 '21+1 '21 

Thus, setting 

4=E/(^U 7 ( fc ))(^)-^)) ) 

Tt '711 ' T 



Z=0 



we deduce that for every k G {0, . . . , 2 k 1 }, 

2 k -l 



h = 4+i + ^2 £ j( k ' 1 ) 



1=0 

with 



Ey(k,l) — (/ 7 (t 2Z + \X 7 (fe+1) ) - f-yiT^+i \x' y (k + 1) ) ) ( -B^fc+l) ~ 

V T 2i T 2l + 1 / \ T 2l+2 T 2l + 1 

Furthermore, J = / 7 (s,X7)(J3f - 5f ) and I* = #7 - fl7- We then deduce that for 
every 7 > 0, 

if-12 fe -l 

H?-iq = f 7 (8 J Xr ) (B* r -£f)-E E £ t(M)- (17) 

fc=0 «=0 



Let us focus on the second part of the right-hand member of (|17p . Owing to the assump- 
tions on / 7 , we have 



MMI < crty ((rg^ -r 2 ( f +1) ) + |x 7 (fc+1) -^VdO 

\ T 2l+1 T 2l J 



n Jk+l) — n (k+l) 



with 4>k,i = 1 + \X\ +1) \ r + \X\ k+1) \ r . Then, 



(fe+i)i ~ l (fc+i) 1 

r 2l T 2Z+1 



V T 2i+1 T 2i + 1 / 

(18) 

for every G Let us now focus on the control of (j>i^ and of |X 7 (fc+1) — X 7 fe+1) |. 

T 2Z+1 T 2l+1 

First, note that if 4+1^ = £ then, = 4+1 = t an d it follows that in this case, 

e y (k,l) = 0. In the sequel, we set At '■= {(&)0>4+l < t}' By Lemma (TJ for every 
(k,l)eA t , 

M<c T [i + \xi\ r + {\\r'\\^y) . 

Second, using the notation n^{k, I) for the discretization index related to r, (i.e. n 7 (/c, /) = 
7 -1 t/^), we can write 



n 7 (2Z+l,fc+l)-l 



x\ k+1) - x\ k+1) = Y, + ^Vd - z 

1 I I TV, 1 TV, . ■ -, 



7 

(fc+1)- 



21+1 '21 ,„,,.., '2i+l '2; 

•u=ri 7 (2f,K+l) 

Then, using that b is a sublinear function, we deduce from Lemma Q] that there exists 
Ct > that does not depend on 7, I and k such that 

n 7 (2/+l,fc+l)-l 

£ 7 \b(x^)\ <c T (i + \xj\ + P 7 iif;; 2i+i ~ 7 ) (4^ - 

u=n 7 (2Z,fc+l) ^ ' 

<C T (l + \jq\ + H^Hf;^) (rg^ - r 2 ( f +1) ). 

Likewise, we have for every (k, I) G A, 

\?7 y~1 I < II yi\\^t-lr (fc+l)^ 

1^ (fc+l) _z (fe+i)l - ll z lle,7 l r 2«+l _T 2Z i • 

r 2i+l T 2i 
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Plugging the previous controls into (|18j) and using that 29 < 1 + 9, we obtain the existence 
of a constant Ct > such that for every 7, I and k, 



MM)I < c T Mxi ll^ 7 llf 7 " 7 ))(((i + \xZ\)(t$$> - T% +1) t 



+ War., — T 9 , ) \\\t> \\ot- 



',7 V'22+2 

where 

ip r (x,r) = 1 + \x\ r + \z\ r . 
Then, using (fT6j) and the fact that 

2 fc -l 2 fc -l 

E(fc+1) _ (fc+1) V"^ ( fc ) _ ( fc ) - + _ 
T 2l+2 T 2l — 2-^i T '+l T l ~ - -' 
Z=0 Z=0 

it follows that 

Mn , ||z,|,#> £ ( (1 + l^fj + W J 

Summing over k and using that X^fc>o 2~^ k+1 ^ e and X]fc>o 2~( fc+1 )( 26 '~ 1 ) are convergent, we 
deduce that 



ESojKMJ <r" ^_<A ( ft J- \ Y"y\\ ( tlA\\ II ^7,1^-7 f 26 ^ I, R H, 



A'-12 fc -l 

:>*-7\ 

fc=o «=o (19) 



E E M*.oi<or*v(^, ii^iif;' 



x ( (1 + - + ||^||f ~ 7 (t - s) 2d ) \\B H 



\e,T- 



A expansion of the right-hand side combined with elementary inequalities and the fact 
that 29 < 1 + 9 yields for every < s < t < T, 

EE MM)I < C T (t_-sr (l + \Xj\ r+1 + (\\Z%p) 1+r ) \\B H \\ e ,T 
fc=o 1=0 ^ ' 

and the result follows from (|17p . □ 

In the following lemma, we make use of Lemma [2] when / 7 (i, x) = o~(x). In this 
particular case, we show below that we can deduce a control of the increments of on 
an interval with random but explicit length r\(uS) (which does not depend on 7). 

Lemma 3. Let 70 be a positive number. Assume that a is a bounded Lipschitz continuous 
function. Then, for every 9 6 (^,H), for every T > 0, there exists Ct > 0, there exists a 
positive random variable 

:= ^[(CtI^^t)- 1 Al^ (20) 
such that a.s for every < s < t < T with t — s<r], for every 7 6 (0, 70) 
\Zl - ZJ\ <(t_- s) e (Hoc + C r (l + \X2\)f] e ) \\B H \\ e>T 
where ||o"||oo = sup^^d ||cr(a;)||. 
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Proof. For every I > 0, set t\ = s + jl and Ni = ||Z 7 |||'^. Owing to the definition of ||.||§' 7 ; 
we have 



Ni +1 <NiV sup 



Z 7 - Z 7 



i<l (.4+1 ~~ Hj 

By Lemma [2] applied with s = U, t = ti + \ and f y (s,x) = a(x) (and r = 0), 



Z 7 - Z 7 

H+l H 



(U+i - Uf 



< ( IklU + C T (t l+1 - UY ( 1 + \U\ + ||Z 7 ||f < ) ) \\e, T . 



By Lemma [T] and the fact that t — > ||Z 7 |||'* is nonincr easing, it follows that 



Z 7 - Z 7 



sup ■ 

i<; (4+i — h 



< (JMloo + C T (Jl + \Xj\)(t l+1 -s) + \\Z^\\f%(t l+1 - s 

Let p be a positive number. If — s < p, we obtain that 

iVz+i < JVj V (a, + P p Nt) 

with 

' kHoo + C r ((l + |X 7 |)p^ ||£ H || e , T and, /L = qr/||B H ||« T . 



R# II 



Let us now set p = rj(u)) where rj{uj) is defined by ()20p . For this choice of p, we have 
Aj < \ • Then, the interval [0, a v /(l — /3 V )] being stable by the function iHa, + P v x, we 
deduce that for every I £ N such that tj +1 — s < r](u), 



i _ 2 



The result follows. 



□ 



3.2 Proof of Proposition [3] 

(i) We first prove (flQ|) for p = 1. Set A n = 5^ — B^,^. Owing to the Taylor formula, 
V(X (n+lh ) = V(X ni ) + 7 (W(X n7 ),6(X n7 )) + {VV(X nj ),a(X nj )A n+1 ) 

+ ^ /^d« 2 j^(£n+l)(^(n+l)7 ~~ ^n 7 )i(^(n+l)7 ~ X n ^)j. 



where C(n+i) £ P^n 7 , X(n+l)y\- Using Assumption (C) and the boundedness of D 2 V and 
a, we obtain 

V{X {n+l)l ) < V{X ni ) + 7 (/3 - a^(X n7 )) + Ai(n + 1) + C( 7 2 F(X„ 7 ) + |A n+1 | 2 ). (21) 
where 

Ai(n+1) = (Vy(A > n7 ), ( r(A > n7 )A n+1 ). 

a 

Set 70 = — For every 7 S (0,70], for every n > 0, we have 
20 

^(X (n+1)7 ) < V(X m )(l - | 7 ) + A x (n + 1) + (y? 7 + C|A n+1 | 2 ). 

12 



Then, iterating the previous inequality yields for every s, t such that s < t, 



7 

C\ t — s _ , r\ i — s , 

V(X t )<V(Xs)(l-- 1 )-r+ £ (1 - (Ai(fc) + /3 T + C|A fc | 2 ) . 

7 

Using that log(l + x) < x for every x > —1, we deduce that 



t — s 



V{X t )<e-^ A {V{X L ) + \Hl-Hl\)+ £ (/3 7 + C|A fc | 2 ), (22) 



fc==+i 

7 



where 



Jo 77 Jo 



with # 7 (s) = (1 - ^)~^- For every G {l,...,d} x set fi' J (s,x) = 

g J (s)(VV)i{x)a ijj (x). Using that sup te[0]T]i7e(0i70 ] \g' y (t)\ < +oo, we check that 0y(-)) 7 e(o, 7o ] 
is a family of Lipschitz continuous functions such that sup 76 ( 0)7o i [<? 7 ]Lip < +oo. Further- 
more, (VV)i and <7jj being respectively Lipschitz continuous and bounded Lipschitz con- 
tinuous functions, we deduce that (/ 7 J ) 7g (o, 7 o] satisfies (JT3|) with r = 1. Applying Lemma 
[2J we obtain that for every 6 6 (|,iT), 



1^7-^71 



(1 + \X]\) + C T (i - s) e (l + \XJ\ 2 + (H^llf^ 7 ) 2 ) 



I R-ff II 



Now, if t - s < T)(u) defined by (|20]>. 

\\Z~<\\lp < (||a|U + C T (1 + |X S |)/) ||£*|| e , T 
but owing to the definition of 77, we have a.s. 

[|B H M|| 9j7 V< C T 
where is a deterministic positive number so that 

(||^|||f T ^) 2 <C' r (||^||^ T + l + |X s | 2 ). 

Thus, 

|F7 - F2| < C T [(1 + - s) e + (t_- s) 29 (l + |X<f + \\B H \\l T )] \\B H \\e, T . 

Using that \ab\ < 2~ 1 (|a| 2 + |6| 2 ) and that 1 + \x\ < CW(x), we have 

(1 + \X2m-s) e \\B H \\ e> T < C{V{Xl){t_-sf e + \\B H \\l T ). 
It follows that there exists Ct > such that for every e > 

\Hl-Hl\<e(t_-s)V{Xl) ^ <^a~ ^"Hl + ''^''^^ ^(H^llg.r+a-^ll^lll.r). 
Now, we choose fj £ G (0, 77) 

C T (fj £ ) 26 -Hl + \\B H L T ) 



< 1. 



e 
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More precisely, we set f} £ = [{C T {l+\\B u \\ ejT )- L e\^ A 77. Thus, we obtain that for every 
< s < t < T such that t — s< fj e , 

\Hl - HJ\ < e(t - s)V(XD + C T {\\B H \\\ T + (t~ s) 29 \\B H \\l T ). (23) 

Then, we can set Eq > in order that there exists 5 > such that 

Va?e[0,l], e~7 <B (l + e x) < 1 - 5x. 

Thus, setting fj := fj eo and plugging the two previous controls in (j22j) . it follows that for 
every k G {!,... , 



A: 7/ 

I ) 



VjXkfj) < V(X (k _ 1)n )(l - Sjj) + C r (l + [|g H [|g |T ) + ^ (/37 + C|A 
An iteration yields for every £{!,..., [-^J}: 



. (fc-i)g 1 



< F(x)(l - <^) fc + C r [|S ff ||| iT £ (1 - 6v) k ~ m 

mfj 

k T 

m=l ; _ (m-l)fj 

~~ 7 

Now, on the one hand, $^j = i(l ~~ Srj) k ~ l < and owing to the definition of f\ (and of rj), 

i 

v '( have fj -1 < C(1 + \\B H Wg 8 ^ 1 ). It follows that there exists a function Pg with polynomial 
growth such that 

Ct\\B h \\% t £(1 - < P(||B*|| fl , T ). 

m=l 

In the sequel of the proofs the index 6 in Pg, which recall the dependance of the polynom 
in 6 is dropped. On the other hand, since (1 — 5rj) k ~ m < 1, we also have 

m ^7 I kfj I I kfj I 

(1 - £ (/3 7 + C|A Z | 2 ) <J> 7 + C\A U \ 2 ) </3kfj + Cj2 I A u | 2 . 

m=l (m— l)g li=l u=l 

~~ 7 

We deduce that for every k E {1,. . . , 

F(A^) < V(x)(l - 5f,) k + P(\\B H \\ e , T ) + CQ 7 (Pf ,0 < t < T), (24) 



where P is a function with polynomial growth and Q7 is defined in the lemma. Applying 

obtain 

V(Xt) < V(x)(l - 5f,)^ + P(\\B H \\ 0iT ) + CQJB*, < t < T). (25) 



this inequality with k = |_?J , we obtain 



Finally, we want to control V(Xt) — V(Xj>). The function VV being sublinear and D 2 V 
being bounded, we deduce from the Taylor formula that for every 1,1/6 M d , 

V(y) < V(x) + C(\x\.\y -x\ + \y - x\ 2 ). 
14 



Applying this inequality with x = Xt and y = Xt and taking advantage of the assump- 
tions on b, we have 

V{X T ) < V{Xt) + C [ 7 (1 + \X Z \ 2 ) + (1 + \Xt\)\B$ - Sf I + \B$ - 5f | 2 ] (26) 
< V(Xt)(1 + Cj) + C(l + \\B H \\% T ), (27) 

where in the second line, we again used the elementary inequality \ab\ < 2 _1 (|a| 2 + \b\ 2 ) 
and the fact that \x\ 2 < CV(x). Combined with (|25|) . the previous inequality yields the 
previous inequality yields: 

V(X T ) < V(x)(l - 5t7) L WJ(1 + C 7 ) + P{\\B H \\e, T ) + CQ^(B t H , < t < T), 

I — I 

where P again denotes a function with polynomial growth. Finally, since (1 - 5fj) l v j < 
e~ ST there exists 70 > and p G (0, 1) such that for every 7 G (0,70), 

(l_^)L|J( 1 + (7 7 ) < p . 

Inequality (fTUI) follows. 



Case p > 1 :For every n,i)6l, \u+v \ p = \u\ p +p\u+ kv\ p l sgn{u+Kv)v with k G [0, 1]. 
Since \u+kv\ p ~ 1 < 2P~ l {\u\ p ~ l + \v\ p ~ l ) , we deduce that \u+v\ p < \u\ p +c p {\v\.\u\ p - 1 + \v\ p ) 
for every u, v G M and p > 1. Then, by the Young inequality, it follows that for every 
e > 0, there exists c £jP > such that \u + v\ p < (1 + e)|^| p + c £ , P \v \ p for every u, v G M and 
p > 1. Applying this inequality, we deduce from the case p = 1 that 

V p (Xl) < + e)V p {x) + C7 (P(||S H || e> i) + CQ 7 (Bf,0 < t < l)) p . 

Since p < 1, we can choose e > such that / o = p p (l + e) < 1. Using again the elementary 
inequality \u + v\ p < 2 p ~ 1 (|n| p + \v\ p ), we obtain 

V p {Xj) < pV p (x) + P p ,e{\\B H \\g tl ) + CQ p (B t H ,0< t < 1) 

where P Pt g is a polynomial function. 

Now, let us focus on (]12p . We only give the main ideas of the proof when p = 1 (the 
extension to p > 1 again follows from the inequality |-u + t>| p < 2 p_1 (|u| p + |u| p )). By 
(|24h . the announced inequality holds taking the supremum of the left-hand side of (I12p 
for every kfj with /c G {1, . . . , |_^J}- Then, for every t G [(k — l)fj, kfj], it remains to 
control (uniformly in k) V(X t ) in terms of V^^^i^). By (|2T]l and (|23|) . we obtain such 
a control for every discretization time between (k — l)fj and kfj. Then, it is enough to 
control uniformly V{Xt) in terms of V{Xt). This can be done similarly as in inequality 

(USD. 

(ii) Let s, t G [0, T] with < s < t < T. We have 

X?-X]= f b(Xl)du + Zl -Z]. 

J s 

First, since |6(x)| < Cy/V(x), 

I f b(X2)du\ <(t-s) sup VV(Xu). 
Js ue[o,T] 
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Second, we focus on the increment of Z 1 . By Lemma [3l for every u,v € [0,T] such that 
v—u<r] (where rj is given by (|20l) ). 

\ZZ-Zl\ < (v-uf (\\a\\ 00 + C T (l+ sup \Xs\)v d )\\B H \\e,T- 

V se [°' T ] / 

Using the concavity of x \— > x d on R + , we have for every s\,82 6 [0,T] being such that 
|S2 - si| < 7, 

\Z] 2 - Zl\ << 2 1 - e ||a|| 00 ( S2 - S i) e ||^|| e ,T 
and we derive that for every u,v G [0,T] with \u — v\ < r/, 

\Z2-ZZ\ <C T (v-u) e ( ||(t[|oo + (1+ sup |X|)r7 fl ) ||£ H ||0, T . 

V «e[Q,T] / 

Now, by the very definition of 77, we have r/H-B^ H^t < 1. Then, since |x| < CV(x), we 
deduce from the first part of this proposition that for every u,v £ [0, T] with \u — v\ < 77: 

\Z2 - Zl\ < C T (v - u) e {V{x) + P(\\B H \\ e>T ) + CQ 7 (Bf ,0 < t < T)). (28) 

where P is a polynomial function and Q 7 is defined in Lemma El 

We want now to make use of the preceding inequality to control Z^ — Z] for every < 
s < t < T. We divide [s,t] in intervals of length lower than 77. More precisely, setting 
Sk = s + k [r]\ , we have 



zl - 2] = zl - zi^ + £ z] k - zi_ x . 



i=l 

Then, we deduce from (1280 that 



\Z< - Z]\ < C T ({t - s^) 6 + L^yV) (V(x) + P(\\B H \\o iT ) + CQ^Bf ,0 < i < T)) 
< Cr ((* - s) e + (t - s)^" 1 ) (U(x) + P(||2? H |k T ) + CQ 7 (Bf , < i < T)). 



Thus, using (j28|) if i — s < 77 or the fact that (t — s)^ 1 < (t — s) e if t — s > 77, we deduce 
that there exists CV > such that for every < s < t < T, 

\ZT t - Z]\ < C T (t - s) e (V(x) + P(\\B H \\e, T ) + CQ 7 (B t H ,0< t < T)). 

The result follows. 



4 Tightness properties 

In the following proposition, we obtain some a.s. tightness results for the sequence (V^'" 1 ^ (uj, da)) n >i. 
Using that the controls established in Proposition [3] are uniform in 7, we also show that 
tightness properties also hold for the set of its limiting measures (W(°°'^(u;,0)) 7 defined 
by 
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Proposition 4. Assume (C). Then, there exists 70 > such that, 

(i) For every 7 G (0,7o] and p > \, a.s., 

1 - 

limsup- rv»(r„ x) ) < C p E[|P p (||^|| M )| ,0 < t < 1)] < +00. 

where C does not depend on 7 and P p is a polynomial function and Q 7 is defined by (jlip . 

(ii) For every 9 G (1/2, H), for every 7 G (0,70], (^"'^(a;, c£a)) n >i almost surely tight 
onC e (R + ,R d ). 

(Hi) For every 9 G (1/2, ff), (W^'T^w, 6 , )) 7 e(o, 70 ] *« in C e (R+,R d ). 

Proof (i) Case p = 1 : We first focus on the sequence (i X^^o 1 ^(-^jfc))iv>i- We set 

Vfc > 0, (5 k B H ) t = B* +t - . 
By Proposition [3] applied with T = 1, we have for every /c > 1 

^) < pT^-i) + + CQ^(k-i)B H )t,0 < t < 1) 

with p G (0, 1). An iteration yields for every A; > 1 

fe-1 

V(Xl) < p k V{x) + ^ p^PiWd^We,!) + CQ 7 ((<5,B H ) t , < f < 1). 



Setting C/i = + CQ^((5iB H ) t ,0 < t < 1) and summing over fe, we obtain 



H |M + CQ 7 ((<^) 
fe=0 v K; fc=0 z=o 



Let us remark that since B H is a C 8 ([0, l^K 9 ) valued Gaussian random variable, the 
norm H-B^H^i has finite moments of every order, which is classical consequence of Fernique 
Lemma. Moreover 

E[Q 7 (Bf ,0 < t < 1)|] < Cj 2H -\ 

Hence 

h(rf) :=E[\P(\\B H \\ eA ) + CQ^B* ,0 < t < 1)|] < +00. (29) 



Since (5iB ) is ergodic and that ^(7) < +00, We have 

N-2 

N 



N-2 

Y / U l ^^E[P(\\B H \\ eA ) + CQ 7 (B^,0<t<l)) a.s. (30) 



1=0 

and it follows that 

N-l 



limsup^Y^ViX^^-^nPil^lle^ + CQ^B^O^tKl)} a.s. (31) 



We want now to use this result to control the a.s. asymptotic behavior of (^ X^fc=o ^(^ 7 fc))«>i- 
By the second point of Proposition [3^i), for every k > 0, 

sup y(^)<C(F(A > ^) + P(||5 fcJ B H || eil ) + C7Q 7 ((^) t ,0<t<l)). 



Je[L*l+i,L , 
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As a consequence, setting N = \j{n — 1)J + 1, we have 

< C( 7 +-)(- + 1) [j= E (^(X fc 7 ) + P(||* fc 5 H ||,,i) + CQ 7 ((6 k B H ) t , < t < 1)) . 

7 V fc=0 / 

Using (pOj) and (j3T|) , the result follows when p = 1 . 

The proof when p > 1 is very similar to the case p = 1 and is left to the reader. 

(ii) By (%), (Pq 71 (w, dx)) is tight on R rf (since V is coercive). Owing to some classical 
tightness results in Holder spaces (see e.g. [27], Theorem 1.4), we deduce that we have 
only to prove that for every T > 0, for every 9 G (1/2, H), for every e > 0, 

1 n 

k=l 

where we recall that 

WT^n (fj(\ \f(t)-f(s)\ 

VT>0, u] 0T (f,S):= sup 



0<s<t<T,0<|i-s|<5 I* ~~ s l 



By Proposition [3] (ii) with 9' G (0,-H"), 



sup B igl < C7 T (t - s ) 9 '- 9 (T/(x) + P(\\B H \\e, T ) + Q 7 (Sf , < t < T) 

0<s<t<T [t - S) 

so that for every s, t G [0, T] such that s < t and i — s < 5, 

sup '^"f/ 1 < C T 5 e '- (V(x) + P(||5 H || fl ,r) + Qj(Bf, <t <T). 

0<s<t<T (t - S) u 

As in (i), this property can be extended to the shifted process: we have for every k > 
sup \ X "k+t- x J k+s \ < CT5 e'-e {v{r ) + Pfll^U, iT ) + Q 7 ((* fc B*) t , < t < T) 

0<s<t<T [t - S) u ' 

By (i) and the ergodic properties of the increments of B H (see (j30|) for similar arguments), 
we deduce that 

hmsup - K e T (Xi (k 1V 8) < C6 e '- e . (32) 

fc=l 



The result follows. 



(in) Let G (1/2, H) and denote by // 7) an element of W (oo ' 7) (u;, 0) and by its 
marginals. By (|29|) and ([31~]) . 

V 7 G(0, 70 ], mJV) < j4"^(7) 
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where p does not depend on 7. Since H > 1/2, sup 7e ( 0)7o ] h(j) < +00. It follows that 
Z//q°°' 7 ^(u;, 9) is a.s. tight in M d (where Uq 00 '^^, 9) stands for the set of initial distributions 

Now, since C does not depend on 7 in (j32l) . we also have for every T > 0, <5 > and e > 
for every 9' > 9: 

V 7 6 (0,70], M (7) (1{^ T M)> £ }) < C5 e '- 9 . 
and the announced result follows again from Theorem 1.4 of [27]. □ 



5 Identification of the weak limits 

5.1 Weak limits of (V {n '^(u, da))„>i 
We have the following result: 

Propositions. Assume (C) and let V^ 00 ^^, da) denote a weak limit of '((V^ n ^(cj, da)) n >\. 
Then, V^ 00 '^ (ui, da) is a.s. an adapted stationary solution of ©. 

Remark 4. In the following proof, we will prove some properties "for every function /, for 
almost every cj" and conclude that "for almost every u, for every function /" the property 
is true. For the sake of completeness, we recall here that such inversions are rigourous 
since we work on Polish spaces (in which the distributions and the weak convergence are 
characterized by some countable family of bounded continuous functions) . 

Proof. In the proof, we denote by {V^ n \uj,da,dj3)) n >i, the sequence of probability mea- 
sures on C e (IR + ,M d ) x C e (R,Ri) with \<0 < H defined by 

1 ™ 

Vin,l) = - E % 7 Vd+. M'^-ib*. M-*5-ib M) (da ' dP)) 

k=l 

where (B^)t^u is the fractional Brownian motion used to build the Euler scheme (|9"j). 
First, let us recall that by Propositions] (ii), {V (jin \uj, da)) n >i is a.s. tight. Thus, we can 
consider a weak limit 'P( 00 ' 7 )(u;, da). Second, one checks that {V^ n,1 '{oJ, da, d/3)) n >i is also 
almost surely tight since each of its margins have this property. Indeed, for the first margin, 
it is again (ii) of Proposition [H For the second margin, we use that (B^^r is ergodic 
under the transformation T 7 : C e (R, R d ) -> C e (R, R q ) defined by (T 7 (w))t = w(7+t)-w(7) 
(see e.g. [22J). In particular, 

1 ™ 

->5 r h rh (dfl) (33) 

k=i 

is converging almost surely to the distribution of {B^) t< ^ (on C e (R, R d )). Hence, the 
sequence almost surely tight (and thus relatively compact). Then, 

if T > ^°° ,n, \u3, da) is the limit of a subsequence of (V^""'^ (uj, da)) n >\, maybe with the help of 
a second extraction, it follows that a.s., there exists a subsequence (nk(u}))k>o such that 

V {rik '^(u,da) k ^ + °°) V^\u,da) and V {nk ^{u,da,dp) " fc ^ +0 °) p<P°>i)(u, da, dp) 

(34) 

where the first margin of V^^'ioj, da, d/3) is obviously 'P'- 00 ' 7 *' (u, da) and the second one 
is a.s. the distribution of (B^)t^R (thanks to ([33]) ) . Let us also denote by {X[°° ,B^) 
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the coordinate process on C 9 (IR + ,IR d ) x C^QR,!^) endowed with the probability p(°°>T). 
For {a, 13) G (7 e (IR + ,]R d ) x C e (R + , R q ) we consider the following function 

$?(a, p) t := a + / H&(a, 0).)da + f 0).)dp.. (35) 



Please remark that <3? 7 is slightly different from <3? 7 in the way it handles the initial con- 
dition but 

= $>,/?) 

for every a such that cto = a. For t, K > let us denote by the function defined 
on C ,e (M + ,M d ) x by F tjK {a,P) = sup < s < 4 |a £ - &{a',/3+) s \ A K where 0+ = 

(P(t)) t > . The function F t>K is bounded continuous on C (R+,R d ) x C e {l 

Then, 



1 

E(*U(* ( °°"\ = Jim - £ ^(X ( V 1)7+ , Bg_ 1)7+ . -Bg_ 1)7 ). 

' fc=l 

By definition of the Euler scheme (even though it is shifted), we have for every k > 1, 
F i,^(^(Vi) 7 +.'- B (Li)7+. - B (k-i)J = almost surely, and 

^(00,7) _ $1 (xt 00 ' 7 ^, B H ) 



almost surely, which ensures that the pair is a solution of (EI). 

The stationarity of _X"(°°' 7 ) follows from the construction. Actually, using the convergence 
we have for every bounded continuous functional F : C (R + ,R ) — > R, 



1 n 

fe=l 

and owing to a change of variable, it is obvious that for every t G 7N, 



1 " 

fc=i 



It follows that for every t G 7N, for every F, 

E[F{X^ rf) )} = E[F(X(°°^)]. 
This property implies that X^ 00 ^ is stationary. 

We now focus on the adaptation of X^ 00 ' 7 **. In this step, we need to introduce, for a 
subset D of M that contains 0, the Polish space W# 5(F)) that denotes the completion of 
C^(D,R q ) (space of C°°-functions / : D ->• with compact support and /(0) = 0) for 
the norm 

imi \f(t)-f(s)\ 

This space is convenient to obtain some Feller properties for the conditional distribution 
of the fractional Brownian motion given its past. More precisely, by Lemmas 4.1 to 4.3 of 
[33], the paths of B H belong a.s. to W fl) j(R) when <9 G (1/2, H) and + <5 G (if, 1) and, in 
this case, for every T > 0, 

Vt{u, .) := £((Bf ' T ) t <o|(#f )t<o = (ut)t<o) 
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is a Feller transition on W^j(M ) where B^' T = B^ T — B^ . Let us now prove that X^ 00 '^ 
is adapted. It is enough to show that for every < t < T, for every bounded continuous 
functionals / : C e ([0,t],R d ) -> R, g : W 9 , 5 ((-oo, T]) -> M and h : W^R-) -»■ R 

E[/(Xj°°"\ fl € [0,i])^(Bf ,s < T)h(B? +t ,s< 0)] 

= E[/(X(°°^, S € [0,t])^(Bf , a < 0)/»(B£ tJ * < 0)] 

where ^((uj s ) s < t ) = E[g(Bf,s < T)\(B^) s < t = (cj s ) s < \. Owing to the Feller property 
on Vt(w, .), we easily obtain that ip 9 is continuous on Wg t s(M.-). 

Then, using the ergodicity of the increments of B H , we can show as in the beginning of 
the proof that (P ( "' 7) (w))„>i is tight on C (R+,M. d ) x W e ,s( R )- Thus > there exists a - s - a 
sequence (n k ) such that 

n k 

E[f(X^\s < t)g{Bf, s > T)h(Bf, s<t)}= lim — V H k ^J k 

fc->+oo n k ^— d 

n k 

K[f(X^\s< t)^(B",s< t)h(B?,s< t)\ = lim — J] H k ^E[J k \^ {k _ 1)+t ] 

K k=l 

with F t = a(B*,s < t), H k = f(X7 (k+s ,s < t)h(B* k _ 1)+a+t - B^ {k _ iy s < 0), and 
J k = g(B^ k _ 1 ^ +s — B^ k _ 1 y s <T). This implies that it is now enough to prove that 

1 U n->+oo 

- ^2 H k-\ {Jk ~ E[J k |J r 7(A ._ 1 ) +t ]) °°> a.s. 

n k=l 

This point follows from a decomposition of the above sum in martingale increments and 
from classical martingale arguments (see proof of Proposition 6 of [3] for a similar proof) . 

□ 



5.2 Identification of limits when 7 — > + 

In this part we fix a //-fractional Brownian motion B H on C e (M, R 9 ) and we consider a pair 
(X°°^,B H ) on C e (R + ,M d ) x C (R+,R q ) such that for each 7 > the joint distribution 
is given by Proposition [SJ 

Proposition 6. Let (7^) be a sequence converging to such that the distributions of 
(X°°^ k ,B H ) are converging weakly on C e (R + ,R d ) x C (R,R g ) to (X°°,B H ). Then the 
pair (X°°,B H ) is a stationary adapted solution to (P) in the sense of Definition^ 

Proof. Let us first introduce 

$(a,P) t :=a + [ b($(a, p) s )ds + [ a{$(a, P) s )d/3 S , 
Jo Jo 

and remark that <&(«, /?) = ^(a, /3), if ao = a - We want to show that 

X°° = $(X°°,B H ) (37) 

almost surely so that (X°°, B H ) is a solution to ([1]). Let us rewrite the equation with the 
help of two continuous operators on C ,e (M + ,R ci ) x C (R + ,R<?) : 

« r (a,/3) t = I b(a s )ds+ [ a(a s )dp s , 
Jo Jo 
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and 

A(a) t = a t - a . 

Then equation (f3"T|) is equivalent to 

AiX 00 ) = ^{X°°,B H ). (38) 
Let us also consider the discretization of ^ 

*T(a,0) t = / 6(a. )ds+ [ a(a s W s . 
Jo Jo 

Obviously © can be rewritten 

A(X°°^) = *T(X°°'T, (39) 

Lemma 4. Let (7fc)fc>i fre a sequence converging to such that (A"°°' 7fc , B H )k>i converges 
weakly on C e (R+,R d ) x C e (R,R<?) to (X°°,B H ). Then W°(X°°>v<, B H ) converges weakly 
on C e (R+,M. d ) to y(X°°,B H ). 

Proof. Let (a,/?) G C ,fl (M + ,M (i ) x C ,e (M + ,R 9 ). A classical result concerning the discretiza- 
tion of Young integrals shows that 

\V{a,P) t -W(a,P)t\ < ||a|| e , t || J 8|| fl , t 7 2fl - 1 t. 

See for instance 0], Proposition 31 or [30]. Hence for T > 0, 

||*(a,/3) - W{a,P)\\e,T < H^tWPW^t^T 1 - 6 . (40) 

Let F be any bounded -fT-Lipschitz functional on C (R+, [0, T]), 

|E(F(^(X°°'^,^)) -E(F(*(X°°, J B H ))| -»• (41) 

as — )• oo. Then 

|E(F(^(^°°' 7fc ^ H )) -E(F(^ (X°°>i\B H ))\ < KEQ\X oo ''»'\\ eT \\B B \\ eT )T 1 - >$- 1 , 

(42) 

and using Proposition [3](ii) the left hand side of (|42p is converging to as — > oo. 
Combining ()4ip and this last fact, we get the desired convergence in distribution. □ 

Let us start with 

A(X°°> 7fc ) = ^(X°°^, B H ), (43) 

and let k — > oo. By Lemma HI the right hand side of (|43p converges to ^>(X°°, B H ) and 
the left hand side to A(X°°), which, in turn, has the same distribution as ^(X°°,B H ). 



Now, let us prove that X°° is stationary. It is enough to show that E[i ? (X°°)] = 
E[F(X t °^)] for every t > and for every functional F defined by F(a) = IJfcLi fi( a u) 
where fi, ■ ■ ■ , f m denote Lispchitz continuous functions on R rf and t\, . . . , t m belong to 
R+. By Proposition [5j the distribution of X 00 ' 7 is invariant by the time-shift (#fc 7 ) for 
every k G N so that E[F(X°°)} = E[F(X^_ )). The result follows easily by checking that 
for every T > 0, 

E[ sup |JC' 7 -X~' 7 |] ^>0. 

«,ue[0,T],|u— y|<7 

Finally, it remains to show that (X°°,B H ) is adapted. Since (X°°' 7fc ) converges in dis- 
tribution to X°° on C e (R + ,R d ) and since B H belongs to W e ,s (with 9 G (1/2, H) and 
9 + 5 e(H, 1)), (X°°^'k,B H ) converges to (X°°,B H ) for 7^ a subsequence of j k . Then, we 
can let 7 go to in equality ([36]) and the result follows. □ 
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6 Simulations 



In this section, we give an illustration of the application of our procedure for a one- 
dimensional toy equation: 

dX t = -X t dt + (4 + coB(X t ))dB? . 

We propose to compute an estimation of the density of the (marginal) invariant distri- 
bution in this case. We denote it by v$ . By Theorem [Q for every bounded continuous 
function / : R d -)• R, 

lim lim V^\co,f) = ^(f). 

■y—>0n— >+oo 

The first step is to simulate the sequence (B^ k — B5 fc _ 1 \)fc=i- We use the method of 
Wood-Chan (see [29]) which is based on the embedding of the covariance matrix of the 
fractional increments in a symmetric circulant matrix (whose eigenvalues can be computed 
using the Fast Fourier Transform). 

Then, we compute * "Pq™' 7 ^ where is the Gaussian convolution kernel defined by 



K h {x) 



exp(- 



2h> 



Note that K h *V^' y) {x ) 



T >< ^ 1, ' 1 \K} 1 {xq — .)), where, for a 



2h,r — "» ' ' U - o 

measure fi, and a /U-measurable function f, we set //(/) = / /<i)U. In Figure [His depicted 
the approximate density with the following choices of parameters 



n 



10 7 , 



7 = 0.05 h = 0.2, H 



3 
4' 



We choose to compare it with the density of the invariant distribution when H = 1/2. 
Note that in this case, the invariant distribution is (semi)-explicit (as for every ergodic 
one-dimensional diffusion) and is given by 



M(dx) 
Af(R) 



where M{dx) 



(4 + cosx) 2 



cxp 



2u 



(4 + cos u) 2 



du dx. 



We observe that the distribution when H = 3/4 has heavier tails than in the diffusion case. 
Finally, in order to have a rough idea of the rate of convergence, we depict in Figure [2] the 




Figure 1: Approximate density of Vq (continuous line) compared with that of v§ (dotted 
line) 



approximate densities for different values of n keeping the other parameters unchanged. 
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Figure 2: Approximate density of for n = 10 5 (dotted line), n = 10 6 (dash-dotted 
line), n = 10 7 (continuous line) 

Remark 5. As mentioned before, this section is only an illustration. In fact, there are 
(many) numerical open questions. For the estimation of the error, it would be necessary for 
a function / to get some rate of convergence results for V < ^ in \f) — vji(f) (long-time error) 
and for v^^if) — fjf (/) (discretization error) where denotes the initial distribution 
of the stationary Euler scheme with step 7. Note that in the diffusion case, the long time 
error is managed by a CLT with rate ('jn)~ a whereas the discretization error is 0(7) 
(see [28j). Finally, even if the Wood and Chan simulation method is fast and exact, it 
requires a lot of memory because of the Fast Fourier Transform. On Matlab, for instance, 
this implies that we can not take n greater than 2.10 . Thus, it could be interesting to 
study some discretization schemes based on some approximations of the fBm-increments 
simulated, which consumes less memory. 

7 Appendix 

Proof of Proposition [JjLet us show that (X~p.) is a skew-product in the sense of [12] 
as follows. For a fractional Brownian B H motion on M, set for every n 6 Z = 
~^(Vi+i)7 — ^nj- Setting W := (M d ) z ~ , we then introduce the regular conditional probability 
V 1 : W -»• Mi(R d ) defined b}0: 

and denote by V 1 the Feller transition on W defined for every measurable function / : 
VV -»• M by V^fiu) = J Rd f(oj U £o)p-t(uj, dCv) where for u G (R d ) z - and Q e R d , oj U Co = 
(. . . , u) 2 , oj 1 , ljq, Co). Setting <£ 7 (x, Co) = x + 76(3;) + a(x)Co and := £((A n ) n <o), we have 
defined a skew-product (W, "P 7 , R , < I )7 ) with the transition operator Q 7 on M ' x W 
defined by 

q 7 /(x, W )= y /($ 7 (x,^ 7 (^'), 

which describes the dynamics of the Euler scheme. 

Then, thanks to Theorem 1.4.17 of [13], uniqueness of the adapted and stationary discrete 
Euler scheme (X^k) (in distribution) holds, if the skew-product (W, P^-, V , M d , < I )7 ) is 

1 Note that since (A^) nSZ is a stationary sequence, C(AJ\(Al) k < = w) = £(A 7 +1 |(A 7 +fe )fc<o = w) 
for every n£Z. 
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strong Feller and topologically irreducible (in the sense of Definition 1.4.6 and 1.4.7 of 
[13]). 

First, write Co = (Co 1 , . . . ,Co q ) and $ 7 = ($J, ... , Denote by M*(x, Co) the (discrete) 
Malliavin covariance matrix of $ defined by 

d 

\f(x,Co) eM. d xR d aad(i,j) e{l,...,d} 2 , M*(x,Co) :=^d&*7(x,Q)d&$](x,Si). 

k=l 

Thus, M®(x,Co) = (o-cr* )(x) and since a 1 is bounded (and continuous), it follows that 
x — > (det(M*) _1 (x, co) is bounded continuous. Second, the functions D u <&, D W D X & and 
-D^<3? are clearly bounded continuous. Finally, the sequence ((A^) 1 ) has a spectral density 
/ that satisfies (f(x))~ 1 dx < +oo (see e.g. [2] for an explicit expression of /). Thus, 
it follows from Theorem 1.5.9 of [13] that the skew-product is strong Feller. 
For the topogical irreducibility, it is enough to show that for every (x,lo) G M. d x W, for 
every (y,e) G R d x M+, Q(x,uo, B(y,e) x W) > . Since a is invertible, the map $ is 
controllable in the following sense: &(x,Co x ) = y has a (unique) solution Co G M 9 , for every 
x, y G Furthermore, 6 and cr being continuous, for every e > 0, there exists r £ such 
that for every Co G B(Co x ,r £ ), <I>(x,tD) G B(y,e). Thus, 

Q(x,w,S(y,e) x W) > V(io, B(Co x ,r £ )) > 

since V(co, .) is Gaussian with positive variance. This concludes the proof. 
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